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ABSTRACT 

Recent observations on the GRB 000301c afterglow show that three breaks 
appear in the R-band light curve, and in particular the decay slope at late 
times is as steep as —3.0. This unusual afterglow is clearly inconsistent 
with the standard afterglow shock model. Here we propose a non-standard 
model for the unusual R-band afterglow of GRB 000301c. In this model, an 
ultra-relativistic shock in a dense medium ("dirty environment") rapidly evolved 
to the non-relativistic phase in initial 1 day. During such a phase, the shock 
happened to be caught up with by two energetic shells ejected from the central 
engine at two different times, and the shock was refreshed, leading to two 
flattenings of the light curve. After each interaction between the shock and shell, 
the afterglow decayed as oc t~ 3 if the electron distribution index of the shocked 
medium, p 3.4, derived from the optical spectrum. Therefore, this model 
can provide an excellent explanation for the flattening and steepening features 
of the GRB 000301c optical afterglow light curve. We further point out that 
the energy injection shells ejected from the central engine at later times may 
be material shells (e.g., in the massive star progenitor models related to black 
holes) or radiation shells (e.g., in the millisecond pulsar progenitor models). 
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1. Introduction 

The standard model of gamma-ray burst (GRB) afterglows assumes that a relativistic 
fireball is decelerating due to interaction with the surrounding medium (for a review see 
Piran 1999). During such a deceleration, a relativistic forward shock forms and then 
produces an afterglow by synchrotron radiation and/or inverse Compton scattering. The 
simplest case of this model is that the surrounding medium is a homogeneous interstellar 
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one with typical density of ~ 1 cm -3 . In this case, an optical afterglow light curve (e.g., 
GRB 970508) can be well fitted by a single power law until several months. However, this 
property, as we will see below, is clearly inconsistent with the peculiar optical afterglow of 
GRB 000301c. 

GRB 000301c was independently detected by the All-Sky Monitor on the Rossi X-Ray 
Timing Explorer and by Ulysses and NEAR of the current Interplanetary Network on 2000 
March 1.4108. The burst itself had a single peak lasting approximately 10 seconds (Smith, 
Hurley & Cline 2000). Its R-band afterglow on March 2.906 UT was first detected by 
UPSO (Masetti et al. 2000). This burst's redshift was measured asz = 2.0335 ± 0.0003 
(Castro et al. 2000) by identifying weak metal lines in the afterglow's optical spectrum. 
According to the published papers (Rhoads & Fruchter 2000; Masetti et al. 2000; Sagar et 
al. 2000), we can see the following features of the optical afterglow: the R-band afterglow 
light curve in about 4 days after the burst was fitted approximately by a power-law with 
an index of ct\ = —0.82 ± 0.20, and in later one day steepened based on another power 
law with an index of a 2 — —3.0 ± 0.10. However, during a period between the fifth and 
seventh days after the burst, the light curve flattened with the third temporal decay index 
of «3 = —0.53 ± 0.50, and subsequently till March 14.60 UT (the lastly observed time), 
steepened again based on a decay index a 4 similar to a<i- In addition, the V-band and 
B-band afterglow could fade down almost simultaneously with the R-band one. 

A successful scenario must explain the flattening and steepening features of the optical 
afterglow light curve of GRB 000301c. To our knowledge, four mechanisms have been 
proposed to account for steepening. First, as the emission comes from slow-cooling electrons 
to fast-cooling electrons accelerated behind a relativistic shock in a homogeneous medium, 
its decay index steepens by a factor of 0.25 (Sari, Piran & Narayan 1998), which is clearly 
inconsistent with the observational result. Second, as analyzed by many authors (Vietri 
1997; Dai & Lu 1998a; Meszaros, Rees & Wijers 1998; Panaitescu, Meszaros & Rees 1998; 
Chevalier & Li 1999, 2000), the afterglow from a relativistic shock in the wind medium must 
decay more rapidly than in the interstellar medium (ISM). For an adiabatic relativistic 
shock in the wind case, an electron distribution index of p ~ 4.3 is required by a large 
decay index of the late-time afterglow of GRB 000301c, «2 ~ «4 ~ —3.0. This would lead 
to a spectral index of (3 ~ —1-7, which is steeper than the observed one, /3 bs = —1.1 ± 0.1, 
derived from the spectrum taken on 2000 March 3.47 UT by Feng, Wang & Wheeler (2000) 
and on March 14.61 UT by Sagar et al. (2000), respectively. Third, the steepening of a 
late-time optical afterglow light curve may be caused by lateral spreading of a jet (Rhoads 
1999; Sari, Piran & Halpern 1999). A difficulty for this mechanism is that the degree 
of steepening found by numerical studies (e.g., Panaitescu & Meszaros 1999; Moderski, 
Sikora & Bulik 2000; Huang, Dai & Lu 2000; Wei & Lu 2000) when two effects such as the 
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equal-time surface and detailed dynamics of the jet are considered is much weaker than the 
one analytically predicted. Finally, we recently suggested that the evolution of a relativistic 
shock in a dense medium to the non-relativistic phase should lead to steepening of an 
afterglow light curve (Dai & Lu 1999). We found that this model is quite consistent with 
the observations on the GRB 990123 afterglow if the medium density is about 10 6 cm -3 . 
Furthermore, as shown analytically and numerically by Dai & Lu (2000a) and Wang, Dai 
& Lu (2000), this model can also well fit all the GRB 980519 afterglow data. 

Energy injection from the GRB central engine to its postburst shock has been widely 
argued to be a plausible scenario causing flattening of an afterglow light curve. This 
scenario can be realized by two different mechanisms: (1) The central engine may eject 
some shells with different Lorentz factors at different times. As the outer shells move 
outward, they begin to interact with the surrouding medium and decelerate, forming a 
forward shock (blast wave). Eventually the slower inner shells catch up with the outer 
shells. The interaction of slow shells with faster shells that have been slowed down implies 
refreshment of the shock, leading to a flattening of the afterglow light curve (Rees & 
Meszaros 1998; Panaitescu et al. 1998; Kumar & Piran 2000; Sari & Meszaros 2000). (2) If 
the GRB central engine is a strongly magnetized millisecond pulsar, its rotational energy 
input to the postburst shock through magnetic dipole radiation also results in a flattening 
of the afterglow light curve (Dai & Lu 1998b, 1998c, 2000a). In this Letter we argue that 
combination of our dense medium model with such an energy injection scenario can provide 
an excellent explanation for the peculiar optical afterglow light curve of GRB 000301c. 



2. The Model 

The dense ("dirty") environment of GRBs has been discussed in the literature. For 
example, collisions of relativistic nucleons with a dense cloud is suggested by Katz (1994) 
to explain the delayed hard photons from GRB 940217. The presence of an iron emission 
line in the X-ray afterglow spectrum of GRB 970508 and GRB 970828 reported by Piro et 
al. (1999) and Yoshida et al. (1999) requires that the ambient medium of these bursts is 
rather dense (Lazzati, Campana & Ghisellini 1999). The steepening of the light curves of 
some optical afterglows (e.g., GRB 990123 and GRB 980519) may be due to the transition 
to the non-relativistic phase. This also requires that the medium density is as high as 
10 6 cm -3 (Dai & Lu 1999, 2000a). The medium with a similar density is invoked by Dermer 
& Bottcher (2000) to resolve the "line-of-death" objection to the GRB synchrotron shock 
model. This work is guided by the optical observations of r\ Carinae (a best-studied massive 
star), whose environment is a dense cloud (Davidson & Humphrey 1997). In addition, 
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dense media may appear in some energy source models, e.g., failed supernovae (Woosley 
1993), hypernovae (Paczynski 1998), supranovae (Vietri & Stella 1998), phase transitions 
of neutron stars to strange stars (Dai & Lu 1998b; Wang et al. 2000), and anisotropic 
supernovae (Wheeler et al. 2000). 

Based on these motivations, we here assume that the surrounding medium is dense 
and perhaps at different times the central engine ejects several shells, some of which are 
relativistic and others are non-relativistic. It is widely believed that a slow shell contains 
more energy than a faster shell which is ejected at an earlier time. Collisions between shells 
with large Lorentz factors give rise to internal shocks which are expected to produce GRBs. 
After then, these merging shells decelerate, leading to a forward shock (blast wave), as 
they sweep up the dense medium. The shock can be thought to be adiabatic unless the 
electron energy fraction of the shocked medium is as large as ~ 1 (Dai, Huang & Lu 1999). 
The Blandford-McKee (1976) self-similar solution gives the Lorentz factor of an adiabatic 
relativistic shock: 7 = 1.0E 5 1 2 /8 n^ 1/8 t d f y /8 [(l + z)/2] 3 / 8 , where E 52 x 10 52 ergs is the total 
isotropic energy, n = n 5 x 10 5 cm -3 is the medium density, t^ay — t/1 day is the observer 
time, and z is the the redshift of the source. This equation implies that, at time 

the shock begins to enter the non-relativistic phase. 

According to equation (1), therefore, we find that if GRB 000301c was in a dense 
medium with density of ~ 10 6 cm~ 3 , its postburst shock would be non-relativistic at a time 
less than 1 day after the burst. If, during such a period, this shock happened to be caught 
up with by an energetic homogeneous shell which had been ejected from the central source 
at some time, then the shock would be refreshed and its energy evolved based on 

£shock = E + (— j— ) f L(t)dt, (2) 



, 1 + Zj J t 

where Eq is the initial energy of the shock, to is the time at which the shell started to inject 
energy, and L(t) is the injection luminosity. Assuming that the shell's velocity, energy and 
width are f s heii, -Esheii an d Ai?, and the shock's velocity is f s hock, we can write the injection 
luminosity approximately as 

L(t) « ^ hell ~^ ShOCk ^hcll « iW (3) 

In writing the second expression, we have assumed i> s h e ii 3> f s hock- The energy of the 
non-relativistic shocked medium can be approximated by 

2tt 2 2 3 

-^shock ~ -^-^shock-^shock^^p 0(1 ^shock-^shock? (4) 
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where .R s hock is the shock's radius, m p is the proton mass and c is the speed of light. 
Assuming that -E s hcii ^ E and the energy which the shock had obtained from the shell is 
much larger than E when t 3> t , combination of equation (2) with equations (3) and (4) 
leads to 

^shock-^shock Vtt- (5) 

Because -R s hock oc v shock t, we easily find 

^shock oc r 2/b . (6) 



In the following we consider only synchrotron radiation from the shock and ignore 
synchrotron self absorption. To analyze the spectrum and light curve, one needs to know 
two crucial frequencies: the synchrotron peak frequency (v m ) and the cooling frequency (z/ c ). 
Unfortunately, these frequencies are dependent on two unknown parameters: the electron 
energy fraction (e e ) and the magnetic energy fraction (e#) of the shocked medium. Even 
so, the optical-band frequency is usually much higher than the v m of a late-time afterglow. 
From equation (6), we find the shock's radius -Rshock ~ (5/3)t> s hock£/(l + z) oc t 3//5 
and the internal field strength B = (47resnm p Vg hock ) 1 / 2 oc t~ 2 ^. The typical 
electron Lorentz factor 7 m [m p /(2m e )]e e (t> s hock/c) 2 oc t -4 / 5 and the synchrotron 
peak frequency v m = ^(eB)/[(l + z)2irm e c] oc t~ 2 . The cooling Lorentz 
factor 7 C = 67rm e c(l + z)/(a T B 2 t) oc B^ 2 ^ 1 with a T being the Thomson 
scattering cross section (Sari, Piran & Narayan 1998) and the cooling frequency 
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[eB)/[(l + z)27im e c] oc B 3 t 2 oc t 4 / 5 . The synchrotron peak flux decays as 



F Vm = (1 + z)N e P Um /(A7rD 2 L ) oc R 3 B oc t 7/5 , where N e = (47v/3)R 3 n is the total number of 
swept-up electrons in the postshock fluid, P Um = m e c 2 aTB /(3e) is the power radiated per 
electron per unit frequency and is the luminosiy distance from the source. According to 
these scaling laws, we further derive the spectrum and light curve of the afterglow 

F = j {vlv m Y^- 1)l2 F Vm OC ^-1)^(12-5^/5 if v < Ue 

V \{v c ly m Y^ t \vlv c Y^F Vm o,v-^H 2 ^ if v>v c , [l) 

where p is the electron distribution index (Dai & Lu 2000a). We note that if p = 3.4, then 
a = (12 - 5p)/5 = -1.0 and (3 = -(p - l)/2 = -1.2 are consistent with the GRB 000301c 
R-band afterglow data both in initial 4 days and during a period between the fifth and 
seventh days after the burst. These data indicate ot\ = —0.82 ± 0.20, 0:3 = —0.53 ± 0.50 
and p o hs = —1.1 ± 0.1, which imply a b s ~ 5/9 bs/6. If the afterglow were radiated by 
fast-cooling electrons in the shocked medium, we would find a = 2(1 — /3), which is clearly 
inconsistent with the observational result. Therefore, the GRB 000301c R-band afterglow 
arose from those slow-cooling electrons in the shocked medium. 
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After the shell had input all its energy to the shock, from equations (2) and (4), the 
shock's velocity evolved as f s hock oc t~ 3 / 5 and thus the spectrum and light curve of the 
afterglow became 

j I/ -(p-l)/2 t (21-15p)/10 if y < Vc 

v K \ v-p/W-W 2 if v > u c , { ) 

where the v c is different from that in equation (7) (Wijers, Rees & Meszaros 1997; Dai & 
Lu 1999, 2000a). We can see from this equation that in the case of p — 3.4, the model's 
decay index a = (21 — 15p)/10 = —3.0 is quite consistent with the observational data of 
the GRB 000301c R-band afterglow both during a period between the fourth and fifth days 
after the burst and at later times, a 2 ~ «4 = —3.0 ±0.1. 



3. Discussion 

We have shown that our dense medium model combined with the energy injection 
scenario can provide a plausible explanation for the unusual optical afterglow of GRB 
000301c. Now we want to discuss the two energy injection mechanisms in some details. In 
the first mechanism, the shells which input their energy into the shock are material shells, 
whose energy source may be the core collapse of massive stars to black holes. As the shock 
decelerates, such a shell eventually catches up with the shock. During such an interaction, 
two additional shocks (a forward shock and a reverse shock) might form and perhaps give 
rise to some observational effects on afterglows. Kumar & Piran (2000) analyzed such 
effects of ultra-relativistic shocks. We will make a detailed analysis for non-relativistic 
shocks and discuss their possible effects on afterglows in a future paper (Dai & Lu 2000b). 
It is interesting to note that a small but discernible variability appears in the GRB 000301c 
optical afterglow light curve in initial 4 days, which could be due to these additional shocks. 

In the second mechanism, the shells which input their energy into the shock may 
arise from a strongly magnetized millisecond pulsar. In this case, t in equation 
(2) and i> s h e ii = c in equation (3). Because the magnetic dipole radiation luminosity 
L(t) oc (1 + t/T)~ 2 , where T is the characteristic spin-down age, L(t) can be thought of 
as a constant for t < T while L(t) decays as oc t~ 2 for t > T. In addition, we define a 
timescale t CI based on t CI = E /L. As found by Dai & Lu (1998b, 1998c), the pulsar energy 
input effect can be neglected for both t < t CT and t > T, but the pulsar inputs its rotational 
energy into the shock at an approximately constant rate for t CT < t < T . This further 
implies AR « c(T — t cr ) in equation (3). Therefore, the energy injection shells in this 
case are radiation shells. It should be pointed out that the GRB source models involving 
strongly magnetized millisecond pulsars include accretion-induced collapses of magnetized 
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white dwarfs (Usov 1992; Blackmail, Yi & Field 1996; Ruderman, Tao, & Kluzniak 2000), 
mergers of two neutron stars if the equation of state for neutron matter is moderately stiff 
to stiff (Kluzniak & Ruderman 1998), phase transitions of neutron stars (Dai & Lu 1998b), 
R-mode-induced explosions in low-mass X-ray binaries (Spruit 1999), and anisotropic 
supernovae (Wheeler et al. 2000). All these models have been proposed to possibly produce 
GRBs with long durations, which are consistent with GRB 000301c. 

There are two flattenings in the optical afterglow light curve of GRB 000301c. To 
interpret this, we require, in the first mechanism, only two energetic material shells which 
caught up with the shock at two different times. In the second mechanism, how did the 
central pulsar input its rotational energy into the shock twice? We envision that a strongly 
magnetized, rapidly rotating newborn neutron star first spun down through magnetic dipole 
radiation and thus input most of its rotational energy to the postburst shock, leading to the 
first flattening of the afterglow light curve. As the neutron star spun down to some extent, 
it would undergo a phase transition to become a more compact star (e.g., a strange star) 
(Dai & Lu 1998b; Wang et al. 2000). This case is somewhat similar to the one discussed by 
Vietri & Stella (1998), where as the neutron star spins down it will collapse to a black hole. 
Assuming that J is the angular momentum at the phase transition, and Ins an d -^mcs are 
the moments of inertia of the neutron star and the more compact star respectively, angular 
momentum conservation implies that the rotational energy of the post-transition star would 
significantly increase, that is, E MC s = </ 2 /(2/mcs) 3> -^ns = ^V(2/ns) because Jmcs ma y 
be much less than J NS . The post-transition star would subsequently spin down through 
magnetic dipole radiation and inject its rotational energy to the postburst shock, possibly 
resulting in the second flattening of the afterglow light curve. 

The infrared (2.1/xm) light curve of the GRB 000301c afterglow is different from the 
R-band, B-band and V-band light curves. The latter present the flattening and steepening 
features but the former shows a well-sampled break in the decay index at t ~ 3.5 days 
after the burst. The early time index at infrared frequency is very shallow (~ —0.1), 
while the late time index is steep (—2.2) (Rhoads & Fruchter 2000). We here argue that 
such infrared and R-band emissions might result from different radiation mechanisms. 
Otherwise, (1) their temporal decay should be almost independent of frequency, which is 
contrary to the observational result; (2) the ratio of their fluxes should approximately be 
(^k'/^r)"^ 1 ^ 2 ~ 3.7 where p = 3.4. However, this ratio derived from the observed data in 
about three days after the burst is about 10. Therefore, the infrared emission was unlikely 
to arise from the non-relativistic shock. In addition, the singly-broken power-law decay of 
the infrared emission is reminiscent of a relativistic jetted shock. However, as analytically 
shown by Rhoads (1999) and Sari et al. (1999), the time index of an afterglow from a 
lateral-spreading jet evolves from a = 3(1 — p)/4 or (2 — 3p)/A to a = —p. The observed 
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early-time index at infrared frequency (a ~ —0.1) leads to p ~ 1.1 or 0.8, implying that 
the late-time index a ~ —1.1 or —0.8, which is clearly inconsistent with the observations. 
Therefore, the infrared emission couldn't arise from a relativistic jetted shock. We note that 
for the GRB 990123 afterglow its infrared light curve is much different from that at R-hand 
frequency (Kulkarni et al. 1999). The observed infrared afterglow emission of GRB 000301c 
might be produced by dust sublimation (Waxman & Draine 2000; Esin & Blandford 2000). 
However, whether or not this possibility is correct needs further analytical and numerical 
studies. 



4. Conclusions 

Many optical afterglows can be well fitted by a single power-law decay, which supports 
the standard relativistic shock model. But, three breaks appear in the R-band afterglow 
light curve of GRB 000301c, and in particular the decay index at late times is as steep as 
—3.0. This unusual afterglow is clearly inconsistent with the standard model. Following Dai 
& Lu (1999, 2000a), we have here proposed a non-standard shock model for the unusual 
R-band afterglow of GRB 000301c. In this model, an ultra-relativistic shock in a dense 
medium ( "dirty environment" ) rapidly evolved to the non-relativistic phase in initial 1 day. 
During such a transition, the shock was by chance caught up with by one energetic shell 
ejected from the central engine at a later time, and the shock was refreshed, leading to the 
first flattening of the light curve. Once the interaction between the shock and shell finished, 
the afterglow started to decay as oc t~ 3 - if the electron distribution index of the shocked 
medium p 3.4 derived from the optical spectrum. One day later, the shock was by chance 
caught up with by another more energetic shell, resulting in the second flattening of the 
light curve. After this interaction, the shock evolved based on the Sedov- Taylor self-similar 
solution without energy injection and the afterglow light curve steepened again. The energy 
injection shells ejected from the central engine at later times may be material shells (e.g., 
in the massive star progenitor models related to black holes) or radiation shells (e.g., in the 
millisecond pulsar progenitor models). 

We would like to thank Drs. Y. F. Huang, Y. P. Wang and D. M. Wei for helpful 
discussions. This work was supported partially by the National Natural Science Foundation 
of China (grants 19825109 and 19773007) and partially by the National Project on 
Fundamental Researches. 



- 9- 



REFERENCES 

Blackmail, E. G., Yi, I., & Field, G. B. 1996, ApJ, 473, L79 
Blandford, R. D., & McKee, C. F. 1976, Phys. Fluids, 19, 1130 
Castro, S. M. et al. 2000, GCNC 605 
Chevalier, R. A., & Li, Z. Y. 1999, ApJ, 520, L29 



Chevalier, R. A., & Li, Z. Y. 2000, ApJ, in press flastro-ph/9908272|) 

Dai, Z. G., Huang, Y. F., & Lu, T. 1999, ApJ, 520, 634 

Dai, Z. C, & Lu, T. 1998a, MNRAS, 298, 87 

Dai, Z. G., & Lu, T. 1998b, Phys. Rev. Lett., 81, 4301 

Dai, Z. G., & Lu, T. 1998c, A&A, 333, L87 

Dai, Z. G., & Lu, T. 1999, ApJ, 519, L155 

Dai, Z. G., & Lu, T. 2000a, ApJ, 537, in press ( |astro-ph/9906TUg| ) 

Dai, Z. G., & Lu, T. 2000b, in preparation 

Davidson, K., & Humphreys, R. M. 1997, ARA&A, 35, 1 

Dermer, C. D., & Bdttcher, M. 2000, ApJL, submitted ( |astro-ph/ 0002301 ) 

Esin, A. A., & Blandford, R. D. 2000, ApJL, in press ( |astro-ph/00D34T5D 

Feng, M., Wang, L., & Wheeler, J. C. 2000, GCNC 607 

Huang, Y. F., Dai, Z. G., & Lu, T. 2000, MNRAS, in press 

Katz, J. I. 1994, ApJ, 432, L27 

Kluzniak, W., & Ruderman, M. A. 1998, ApJ, 505, L113 
Kulkarni, S. R. et al. 1999, Nature, 398, 389 
Kumar, P., & Piran, T. 2000, ApJ, in press (|astro-ph/9906002| ) 
Lazzati, D., Campana, S., & Ghisellini, G. 1999, MNRAS, 304, L31 
Masetti, N. et al. 2000, A&A, submitted (|astro-ph/000418q ) 



-10- 



Meszaros, P., Rees, M. J., & Wijers, R. A. M. J. 1998, ApJ, 499, 301 

Moderski, R., Sikora, M., & Bulik, T. 2000, ApJ, 529, 151 

Paczyriski, B. 1998, ApJ, 494, L45 

Panaitescu, A., & Meszaros, P. 1999, ApJ, 526, 707 

Panaitescu, A., Meszaros, P., & Rees, M. J. 1998, ApJ, 503, 315 

Piran, T. 1999, Phys. Rep., 314, 575 

Piro, L. et al. 1999, ApJ, 514, L73 

Rees, M. J., & Meszaros, P. 1998, ApJ, 496, LI 

Rhoads, J. 1999, ApJ, 525, 737 

Rhoads, J., & Fruchter, A. S. 2000, ApJ, submitted (|astro-ph/0004057|) 

Ruderman, M. A., Tao, L., & Kluzniak, W. 2000, ApJ, submitted ( |astro-ph/ 0003161) 

Sagar, R. et al. 2000, |astro-ph/000"422l 

Sari, R., & Meszaros, P. 2000, |astro-ph/000"3406 

Sari, R., Piran, T., & Halpern, J. P. 1999, ApJ, 519, L17 

Sari, R., Piran, T., & Narayan, R. 1998, ApJ, 497, L17 

Smith, D. A., Hurley, K., & Kline, T. 2000, GCNC 568 

Spruit, H. C. 1999, A&A, 341, LI 

Usov, V. V. 1992, Nature, 357, 452 

Vietri, M. 1997, ApJ, 488, L105 

Vietri, M., & Stella, L. 1998, ApJ, 507, L45 

Wang, X. Y., Dai, Z. G., & Lu, T. 2000, MNRAS, in press ( |astro-ph/ 99 Iggjg ) 
Wang, X. Y., Dai, Z. G., Lu, T., Wei, D. M., & Huang, Y. F. 2000, A&A, in press 
Waxman, E., & Draine, B. T. 2000, ApJ, in press ( |astro-ph/990902q) 
Wei, D. M., & Lu, T. 2000, ApJ, in press 



- 11 - 



Wheeler, J. C, Yi, I., Hoflich, P., & Wang, L. 2000, ApJ, in press 
Wijers, R. A. M. J., Rees, M. J., Meszaros, P. 1997, MNRAS, 288, L51 
Woosley, S. 1993, ApJ, 405, 273 
Yoshida, A. et al. 1999, A&AS, 138, 433 



This preprint was prepared with the AAS IATeX macros v4.0. 



